Ab initio study of helium and hydrogen interactions in a-Fe 



C. J. Ortiz, 1 R. Vila, 1 and J. M. Pruneda 2 'E 

1 Laboratorio Nacional de Fusion por Confinamiento Magnetico, 
CIEMAT, Avenida Complutense, 40, 28040, Madrid, Spain 
2 Centre d' 'Investigation en Nanociencia i Nanotecnologia (CSIC-ICN). Campus de la UAB, E-08193 Bellaterra, Spain 

(Dated: March 5, 2013) 

Density Functional Theory (DFT) calculations show a weak interaction between hydrogen and 
helium in iron, in contrast to previous reports of a strong trapping of hydrogen at helium^. The 
strong preference of He and H to occupy regions with low electronic density (such as vacancies) 
explains this discrepancy, with vacancy-He and vacancy-H binding forces concealing the repulsive 
interaction between He and H. Furthermore, Rate Theory simulations based on our DFT-calculated 
V n He m H p cluster energetics predict, as it is observed in some experiments, that synergetic effects 
could be expected between H and He in iron under irradiation. 
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Hydrogen and helium are responsible for more than 
98% of nuclear matter in today's universe. Understand- 
ing how they interact is critical for fundamental fields of 
research such as modeling the evolution of stars ^ White 
Dwarfs and the interior of Jovian planets j 4 ^; but also in 
more technological applications such as sequential H/He 
implantation in the semiconductor industry^, or to pre- 
dict the agglomeration of defects and H and He that arise 
by nuclear transmutation between high energy neutrons 
and lattice atoms in materials irradiated in fusion con- 
ditions. In this sense, substantial research efforts have 
been devoted to the experimental and theoretical study 
of the interaction of He-defects, 8 22 or H-defects 23,24 in 
iron-based alloys and different ferritic steels which are 
important structural materials for fusion reactors. A few 
experimental investigations evidenced that helium could 
act as trap for hydrogen, hence reducing diffusivity and 
increasing embrittlement in iron alloys, 25 27 although ex- 
periments could not clarify whether the origin of these 
synergetic effects was due to the interaction between H 
and He itself or due to the trapping at damage produced 
during irradiation. A detailed study of interactions be- 
tween hydrogen, helium and displacement damage is thus 
necessary to understand and predict their simultaneous 
evolution. 

To investigate the origin of these correlation effects be- 
tween H and He in a-Fe is the purpose of this paper. We 
provide an atomistic description of an extensive set of 
small clusters of He, H and vacancies (V n He m H p , with 
n,m,p up to 6), presenting first principles calculations to 
understand the interaction between He and H in iron. 
By comparing the energetics of these clusters we argue 
that the direct interaction between He and H is marginal 
(mostly repulsive) and that the ostensible strong attrac- 
tion reported by others 2 is due to the preferential clus- 
tering of helium or hydrogen to iron vacancies and not 
to the direct interaction between helium and hydrogen. 
Additional kinetic simulations based on a Rate Theory 
approach are presented to illustrate the expected effects 
due to the combined presence of H and He on their evo- 



lution under typical irradiation conditions in a fusion en- 
vironment. 

Calculations are performed using the spin-polarized 
generalized gradient approximation (GGA) 28 as imple- 
mented in the SIESTA code^ Norm-conserving pseu- 
dopotentials 30 with cutoff radii of 1.11, 1.21, 0.71, and 
1.06 A are used for Fe 4s, 4p, 3d and 4/, and a cut- 
off radius of 0.45 A for the partial-core correction. H 
and He pseudopotentials are built with cutoff radii of 
0.66 and 0.69 A, respectively. Valence electrons are de- 
scribed with a linear combination of strictly localized nu- 
merical atomic orbitals that vanish outside certain cut- 
off radii (r c , in A). Triple-^ basis was used for Fe 4s 
(r c = 4.84) and 3d (r c = 3.06) orbitals and double-^ 
for Fe 4p (r c = 3.53), H and He Is (r c = 2.63). 31 Ad- 
ditional (single-^) polarization orbitals were considered 
for the H and He electrons. Supercell calculations were 
first performed on ferromagnetic Fe54 boxes to explore 
the energetics of the different defect configurations, and 
then repeated for larger Fei28 cells. A real space mesh of 
0.079 A was used, and at least 4x4x4 k-point grid sam- 
plings for full relaxation of the atomic positions (forces 
smaller than 0.04 eV/A), keeping the volume of the cell 
fixed. 

The zero-point energy (ZPE) is known to be im- 
portant for H defects, 24 and can be computed from 
the vibrational frequencies of the normal modes using 
ZPE=| J2u where the frequencies v are obtained as- 
suming the Fe atoms are fixed at their equilibrium posi- 
tions (Fe mass is much larger than H or He masses), and 
only the light atoms vibrate (Einstein approximation). 
The ZPE of ^H 2 is 0.13 eV, in good agreement with pre- 
vious results ^ 23 ! 32 and comparable to the average ZPE 
for VH n complexes (0.14 eV/H) reported by Tateyama 
and Ohno*2i Similar values are obtained here for inter- 
stitial H in bulk Fe (0.15 eV and 0.19 eV for octahedral 
and tetrahedral sites), and close to a substitutional He 
(0.15 eV). This homogeneity of the ZPE for H in different 
configurations allows to take an average contribution for 
ZPE of ~ 0.15eV per H atom for each complex cluster. 24 
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FIG. 1: Binding energies of interstitial H to a VH X (gray 
squares), VHeH^ (black squares), and VH^Hx (black dia- 
monds) defect clusters. Also shown the binding energies for 
VH X reported by Tateyama et. al. ref.[2j] (light gray cir- 
cles) and the binding energies to a single substitutional He in 
tungsten (empty squares)^. The horizontal axis refers to the 
number of H atoms in the initial cluster, x. 



One approach to quantify the interaction between H 
and He is through their binding energies. The binding 
energy of hydrogen to a defect cluster, V n He m H p , is de- 
fined as the energy gained by trapping an additional H 
into that cluster to form V n He m H p+ i. A positive energy 
indicates an attraction of H to the defect cluster. Based 
on this representation Jiang and collaborators^ suggest 
a strong hydrogen trapping at helium in W (shown in 
figure [1] with empty squares), with up to twelve H atoms 
clustering around one single substitutional He. Trapping 
at ViHei in iron is notably weaker, with a binding en- 
ergy of 0.2-0.4 eV with a maximum of five hydrogens 
around a helium atom. We found that the binding of 
H to a single vacancy site (no He) is also positive, in 
good agreement with the results reported by Tateyama 
and coworkers 24 (light gray circles) that showed that a 
vacancy can accommodate no more than six hydrogen 
atoms. Notice however the decrease of ~0.1 eV with re- 
spect to the H binding energy to a single vacancy (gray 
squares) observed when a helium atom is placed at the 
vacancy (black squares). 

We can similarly plot the binding energy of an inter- 
stitial He to form a V n He m+ iH p cluster (figure [2]). In 
agreement with previous reports, 11 there is a strong self- 
trapping of He atoms in bulk He (circles in the figure) and 
an even larger trapping at vacancy complexes. However, 
the attractive interaction is again decreased (by tenths of 
eV) when hydrogen atoms are present. In other words, 
although interstitial He lowers its energy by approaching 
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FIG. 2: Binding energies of interstitial He to VnHe m H p com- 
plexes. Empty symbols show the binding to interstitials He 
(circles), VHe x (squares) and V2He x (diamonds) reported in 
ref.jlll]. Our results are shown as solid black symbols. Dark 
(light) gray symbols correspond to the binding energies when 
one (two) H atoms sit close to the vacancy position. The hor- 
izontal axis refers to the number of He atoms in the initial 
cluster, x. 



a vacancy site, the presence of H at the vacancy inhibits 
He trapping. Further evidence that He and H do not 
interact strongly is given by the small (even negative) 
binding energies of He to interstitial H (gray circles). In 
what follows, we will argue that the positive binding en- 
ergies for H at V n He m clusters (Fig. [1]) or He at V n H p 
clusters (Fig. [2]) are not due to the presence of He or H 
atoms, but to the vacancy itself. 

The interatomic potential for He-H in iron can be de- 
termined from the dependence of their binding energy 
on the He-H distance. We can do this for a relevant 
subset of the defect configurations studied in this work, 
which is shown in Fig. [3] together with their correspond- 
ing He-H interatomic distances and formation energies. 
The defect formation energies are defined relative to the 
chemical potentials taken from bulk a-Fe, molecular hy- 
drogen, and the isolated He atom. At the right panels, 
(a) to (d), different configurations for interstitials of He 
and H are presented (VoHeiHi clusters). Left panels, (e) 
to (i), show several ViHei Hi and ViHei H2 clusters, and 
illustrate how hydrogen affects the stability of VHe, dis- 
placing He from the substitutional site, and moving itself 
from the preferential tetrahedral occupation of intersti- 
tial H, towards the octahedral site. 

Central panel in Fig.[3fj) shows the interatomic poten- 
tial (binding energies) for He-H pairs computed in bulk 
a-Fe. Filled squares are obtained from the difference be- 
tween the energy of a HeH complex, and the sum of the 
energy of interstitial hellium and interstitial hydrogen. 
On the other hand, empty squares are defined from the 
decrease (per H atom) on the binding energy of one He 
atom to a vacancy due to the presence of hydrogen atoms 
at the vacancy (difference between black, dark gray and 
light gray squares at x=l in Fig.[2j divided by the number 
of H atoms). This rough approximation follows with rea- 
sonable accuracy the interatomic potential for He and H 
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FIG. 3: (color-online) Atomic configurations for He and H interstitials (a-d) ViHeiHi (e-f) and ViHeilHb defects (g-i) in bulk 
Fe. Blue (dark gray) spheres, orange (light gray) and red square correspond to Fe, H, and He atoms. Vacancy sites are 
marked with a small cross. The interatomic distance and defect formation energies are also shown. Central panel (j) shows 
the interaction energy for He-H pairs in vacuum (solid line) as a function of the interatomic distance, compared to the energies 
for defects in bulk iron. Solid squares are obtained from interstitial pairs and empty symbols are obtained from ViHeiHi and 
ViHeiH2 clusters, as highlighted by dashed arrows. 



computed in vacuum (solid line). It shows that the bind- 
ing of He and H in bulk iron is only marginal, and could 
be modeled with the simple two-body potential of iso- 
lated He and H atoms. Although there is a small attrac- 
tive interaction for interstitials in configurations (c) and 
(d), for which the respective He-H distances are 2.21A 
and 2. 86 A the interaction is clearly repulsive for inter- 
atomic distances smaller than ^2A, which is the case for 
all the studied clusters where He and H are trapped in a 
vacancy. 

These ab initio calculations show that, although the 
direct interaction between H and He is weak, they both 
bind strongly to vacancies, which could be at the origin 
of the experimentally observed synergic effects j^ 5 . - — To 
investigate this point, we simulated the evolution of H 
and He in a-Fe under irradiation, in conditions similar 
to those expected in fusion environment. Nuclear trans- 
mutation continously generates H and He that evolve 
in the presence of a large amount of defects created by 
atomic displacements due to collisions between energetic 
neutrons and lattice atoms. Fusion conditions expected 
in a DEMO-HCLL-4000MW reactor configuration were 
considered. A neutron flux of 1.52 xlO 15 /cm 2 /s was 
calculated 33 and transmutation rates of 1.75 xlO 12 and 
1.1 xlO 13 /cm 3 /s were obtained for He and H in Fe using 
the nuclear inventory code ACAB^i and the EAF 2007 
library 35 . The generation rate of atomic displacements 
(Frenkel pairs) corresponding to the neutron spectrum 
was calculated using the methodology presented in Ref. 
H and was found to be about 1.82 xlO 16 /cm 3 /s. To 
simulate the evolution of H, He and defects we expanded 
the Rate Theory model developed by Ortiz et a& to pre- 
dict the diffusion and clustering of He in a-Fe in the 
presence of impurities. The various V n He m H p clusters 



studied from ab initio in the present work were added 
to the model, with their corresponding binding energies 
and reaction rates, calculated similarly to what is done 
in Ref. [9|. It is important to note that the growth of 
large bubbles containing both H and He is not yet imple- 
mented in our model, since their energetics is not known, 
and only the limited number of small V n He m H p clusters 
is included. 

For our purpose, we considered a temperature of 
450° C, which falls in the range of operating temperature 
(350-500°C) expected for DEMO reactor 37 , depending 
on the cooling options. Fig. [4] reports the evolution of 
the mean bubble size as a function of irradiation time 
for the H-He system described above. For comparison, 
the evolution of the mean bubble size calculated for a 
system where only He is generated (no H) is also shown. 
Our simulations clearly predict that a larger mean bub- 
ble size is expected when H and He evolve simultaneously. 
This result follows the trend of the experimental measure- 
ments obtained in FeCr alloys by Tanaka et. a/., — who 
observed larger bubbles when H is present. It is impor- 
tant to note here that only the mean size of He bubbles is 
shown in Fig. 3J The contribution of the small V n He m H p 
clusters to the mean bubble size is minor. Hence, the 
larger mean size of bubbles that is predicted by our sim- 
ulations when H is present is not due to large bubbles 
containing H but to a synergic effect between H and He. 
Remarkably, a significant effect is predicted though only 
few V n He m H p configurations are included in the kinetic 
model. 

Our simulations confirm that significant synergic ef- 
fects between H and He could be expected in some cases, 
though their direct interaction might be weak, the ap- 
parent strong positive interaction reported in the litera- 
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FIG. 4: Evolution of the mean bubble size in irradiated Fe for 
a system where only He is present (red dots) and for a system 
where H and He evolve simultaneously (black squares). 



ture is due to the much larger V-H and V-He attraction 



that conceals the repulsive interaction between He and 
H. They both tend to be trapped by voids, and this, 
rather than the He-H interaction, could explain the re- 
duction of hydrogen permeation, and the larger cavity 
sizes in ferritic alloys simultaneously irradiated with H 
and He^~— This has important implications for model- 
ing damage evolution in materials for fusion, that until 
now have mainly considered the role of isolated He or 
H aggregates, and not their coexistence. To understand 
in which conditions synergetic effects could be expected, 
further work will be needed on the energetics of larger 
clusters or bubbles containing both He and H in steels, 
and on the effect of solutes and grain boundaries where 
bubbles tend to nucleate. 
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